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Abstract 

A weak-field solution of Einstein's equations is constructed. It is generated by a circular cos- 
mic string externally supported against collapse. The solution exhibits a conical singularity, 
and the corresponding deficit angle is the same as for a straight string of the same linear 
energy density. This confirms the deficit-angle assumption made in the Frolov-Israel-Unruh 
derivation of the metric describing a string loop at a moment of time symmetry. 
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1. Introduction 



In a gauge theory with gauge group G spontaneously broken via the Higgs mechanism to 
a subgroup H, the Higgs field must lie on the manifold of degenerate vacua G/H on the 
circle at spatial infinity Thus if tt\{G/H) is non trivial, string-like configurations with finite 
energy per unit length will exist for which the Higgs field at spatial infinity winds non- 
trivially around G/H. Such so-called cosmic strings have received much attention in recent 
years, especially as providing possible mechanisms for the formation of large-scale structure 
in the universe ||J. Accretion wakes have been proposed to form behind infinitely long 
relativistic strings, which may help explain the observed large-scale filaments. Moreover, 
Vilenkin has suggested that long-lived closed loops of string might act as seeds for galaxy 
and cluster formation. 

The derivation of the gravitational field outside a cosmic string, in the weak-field 
approximation, was carried out by Vilenkin ||, who investigated the case of an infinitely 
thin straight string. The main features of this solution are that the spacetime exterior to 
the string is flat, and that all along the string, the spacetime has a conical singularity with 
deficit angle 5(f) = 87rG/i, where /i is the linear energy density. (The solution to the exact 
Einstein equations have since been derived |3J.) 

Closed string loops have been studied by Frolov, Israel, and Unruh (FIU) || . Using 
the initial- value formulation ||, they constructed a family of momentarily stationary circular 
loops, which are considered as thin loops of string either at the time of formation or at the 
turning point between expansion and collapse. One of the characteristics of the FIU method, 
which has the merit of taking into account the exact non-linear field equations, is that it 
imposes a priori that all points on the circular string be conical singularities of spacetime 
with deficit angle 5<p = 8irGfj,. This is a reasonable assumption because of the observation 
|| that an infinitely thin circular string, when viewed from arbitrarily close to the core, 
is indistinguishable from an infinitely thin straight string, and that all the points on the 
circular string may thus be assumed to be conical singularities with the same angular deficit 
as for a straight string of equal linear energy density. 

The question we ask is the following: instead of making the above assumption, can 
one deduce from the field equations, at least in the weak field limit, the fact that all the 
points on the circular string are conical singularities? As a first step in this direction one 
may construct, without any special hypothesis on conical singularities, a solution produced 
by a circular source generated by a stress-energy tensor obtained by adapting the method 
used by Vilenkin || for a straight string. This solution could then be examined for possible 
conical singularities and the corresponding angular deficits, if any, could be calculated. We 
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carry out such an investigation in this paper. 

An important difference between our analysis and that of FIU is that we seek a solu- 
tion which is stationary, whereas FIU consider a circular string loop which is momentarily at 
rest. Conceivably, a non-stationary solution at a moment of time symmetry may differ from 
a stationary solution by the presence of gravitational waves, but FIU considered precisely 
a class of particular solutions devoid of such free gravitational radiation. Therefore, the 
comparison between the two cases is physically reasonable. Of course an isolated stationary 
circular string is unphysical, since it would tend to collapse. To overcome this problem, 
radial stresses, the values of which are determined by stress-energy conservation conditions, 
are introduced to support the string. It will be seen that these stresses do not contribute to 
the value of the angular deficit. 

In section two, we establish the appropriate form of the stress-energy tensor for this 
problem. We derive and solve the weak-field Einstein equations in sections three and four. 
In section five, we demonstrate that the solution does indeed exhibit conical singularities 
with the same angular deficit as for a straight string with the same linear energy density, 
thus fully supporting and agreeing with the FIU hypothesis. 

In this paper, we use units in which h — c — 1, take the metric to have signature 
(— , +, +, +) and adopt the geometrical conventions of [0]. 

2. The Stress-Energy Tensor and the Metric 

In this section, we establish the stress-energy tensor for a loop of cosmic string arising 
from a spontaneously broken gauge theory. In the case of an infinitely long straight string, 
all components of the stress-energy tensor are localised on the core of the string, that is the 
region of spacetime where the Higgs field unwinds and thus does not lie on the manifold of 
degenerate vacua. The radius of this region is of the order of the Compton wavelength of the 
Higgs field fH, generically a microscopic distance. Thus, when considering the gravitational 
effects of the string, it is physically reasonable to make the thin-string approximation where 
the stress-energy tensor is localised on an infinitely thin line. 

A circular source produces an axially symmetric gravitational field, and the most 
general stationary metric of this type may be cast in the form 

ds 2 = -e 2v dt 2 + e 2C -~ 2v r 2 d<\> 2 + e 2r <~ 2v (dr 2 + dz 2 ) , (2.1) 

where the three functions v, (, and rj depend only on r and z, and x a = (t, 0, r, z) denotes 
cylindrical coordinates. 
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rp V 



(2.2) 



The form of the Einstein tensor for this metric, together with the Einstein field 
equations = —kT^ (k = 8ttG) implies that the most general stress-energy tensor 
compatible with ( |2.1| ) is 

/ a \ 
7 
A e 
V e p J 

where a, j3, 7, A, and e are solely functions of r and z. 

In passing to the limit of an infinitely thin loop of string with radius a, lying in the 
x-y plane and centred at the origin, we obtain (with the exception of T r r , dealt with below) 
the stress-energy tensor as 

V = f ll v 6(r-a)6(z) , (2.3) 

where is the cross-sectional integral of ' . Therefore this limit yields, in an obvious 
notation: 

a = a8(r-a)8(z) (2.4) 

f3 = P5(r-a)5(z) (2.5) 

7 = *f6(r-a)6(z) (2.6) 

e = e5(r-a)5(z) . (2.7) 

However, such a localisation does not occur for the radial stress T r r = A required to support 
the ring. Given that the ring is externally supported, this radial stress is not confined to 
the core. (In what follows, we choose to support the ring from infinity.) 

Following Vilenkin, we can show that (3 = e = 0. Indeed, taking an arbitrary string 
cross-section M. = {((p,r,z) : (ft = const} and using stress-energy conservation iT^\ v = 0), 
we have 

= / T li \x x drdz . (2.8) 
Jm 

In the weak-field limit, the components of the stress-energy tensor ( f2.2| ) are taken to be of 
the same order of magnitude as the metric functions z/, (, and 7], which are assumed to be 
small compared to unity. Thus to lowest order, the only Christoffel symbols for the metric 
Q2.1| ) that contibute to the covariant derivative of the stress-energy tensor are T r ^ = —r 
and = r^0 r = 1/r. Taking fi = z and noting that Tj = = 0, we may integrate by 
parts to get 

= T/ 5 r A + f z z 5 X Z , (2.9) 

(no sum over r and z). There is no boundary contribution since T z x is localised on the 
string core. We emphasise that, in contrast with the case of the straight string, a similar 
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argument would fail to show that T r r = because T r r = A is non-zero outside the core. We 
obtain the desired result (e = ft = 0) by taking A = r, z in (|2.9|) . The last undetermined 
function, A, in the stress-energy tensor can be found in terms of 7 from the stress-energy 
conservation equations, and will be calculated in the next section. 

At this stage, we have deduced that the stress-energy tensor for a thin circular matter 
source supported by external radial stresses is given by ( |2.2|) with ft — e — 0, and a, A and 
7 as above. We rewrite a = —fi and 7 = k, where /i and k are the linear energy density 
and the azimuthal stress per unit length, respectively. Henceforth, we make the thin-string 
approximation and use this form of the stress-energy tensor. 

To specialise to the case of a cosmic string in a spontaneously broken gauge theory, 
we would apply the equation of state || for string matter k = — /i. (This equation of state is 
dictated by Lorentz invariance in the straight-string case; for a circular string, the azimuthal 
stress plays the same role as that of the longitudinal stress T z z for a straight string.) 
We will, however, keep fi and k as two independent parameters. This has the advantage 
that we will then be able to compare and contrast ordinary matter, given by k positive and 
small compared to ji, with string matter, given by k = 



3. The Field Equations and Stress-Energy Conserva- 
tion 

In the weak-field approximation, the non-zero components of the Einstein tensor || for the 
metric ( |2.1|) are, upon retaining only first-order terms in z/, ( and 77: 

Gt = 2V 2 z/ - V 2 C - -C - V 2 r? (3.1) 

r 

G/ = -V 2 r? (3.2) 

Gr = ~Czz Vr (3-3) 

r 

Gz = -trr + ~(Vr- 2 Q (3.4) 

r 

G r * = Czr + ~ (Cz ~ Vz) , (3-5) 

r 

where V 2 = <9 2 + <9 2 , and V 2 = <9 2 + -d r + <9 2 is the flat-space Laplacian operator. (We 
stress that by "weak-field approximation," we mean retaining only first-order terms in u, r\ 
and £.) After a short calculation, the field equations reduce to 

V 2 u = 4vrG(-a + 7 + A) (3.6) 
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V 2 r/ = 87rG 7 (3.7) 
V 2 C + -Cr = 8ttGA (3.8) 



r 



r] r = r( rr + 2( r (3.9) 
Vz = r( zr + ( z . (3.10) 



We integrate Q) and to get 

V = |;(0 + % , (3.11) 

where r/o is an arbitrary constant to be determined later. 

The conservation of energy and momentum [or equivalently the compatibility condi- 
tions for (3.7)-( 3~ToD l places constraints on the allowable components of T^ u , which in the 
weak-field limit read 

7 = |-(rA) . (3.12) 

The solution of (p. 12 ) for A, with the azimuthal stress 7 given by (|2.6| ) and corresponding 
to supporting the string from infinity, is 

A(r,z) = - e(r-a) 6(z) , (3.13) 
r 

in which denotes the Heaviside step function. Furthermore, ( |3.11| ) and ( |3.12| ) imply that 
( p.7| ) and ( |3.8| ) are equivalent. We may therefore discard (|3.8|), retain ( |3.7| ) to find 77, and 
use ( |3.11D to express ( in terms of 77. 

Combining all the results, the field equations are equivalent to solving the following 
set of equations for v, 77 and (: 

V 2 z/ = ATrG[(fi + k)5(r-a) + -Q(r - a)]S(z) (3.14) 

V 2 77 = MGk5{r -a)5(z) (3.15) 

d 

_(rC) = 77 — 770 . (3.16) 

These equations involve the dimensionless quantities Gfi and Gk. For a cosmic string, 
k = —fi, whereas for ordinary matter, k « fi. The largest value of fi in a physically 
relevant theory occurs for GUT strings, for which G^i ~ 10~ 6 0, validating the weak-field 
approximation as an expansion of z/, (, and 77 in powers of Gfi. 
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4. Solution of the Field Equations 



We note that ( |3.14j ) is Poisson's equation for the Newtonian potential v. The solution is 
most easily found in toroidal coordinates (4>,a,ip). The latter are related to cylindrical 
coordinates ((f), r, z) by 

z = aN~ 2 sinip r = aN~ 2 sha (0 < a < oo , \ip |< vr) , (4.1) 

where N 2 = N 2 (a, ip) = cha — cosip. The surfaces o = a = const, are tori whose generating 
circles have radii acosechcxo and acothcio. In particular, the ring r = a, z = is now given 
by u = oo. (See Fig. 1 in ||.) In toroidal coordinates, the metric fl2.1|) reads 

ds * = _ e 2u dt 2 + a 2 N ^ ^ e 2 C -2, # 2 + fl 2 ^-4 £ 2 V -2v ^2 + #2) ? ^ 

where u, (, and rj depend on a and ip. 

The potential v can be split into two parts, v = v\ + z/2, where z/i and satisfy 
separate Poisson equations: 

V 2 v l = 47rG(iJ + k)5(r-a)5(z) (4.3) 

V 2 i/ 2 = 4vrG-e(r-a)5(z) . (4.4) 
r 

The function v\ is the gravitational potential produced by a circular ring of matter. This 
problem is known and the corresponding solution is 



v x = -2 3 ' 2 G([i + k)N(*,4)K(th{±a))/<h(^o) , (4.5) 

where K denotes the Complete Elliptic Integral of the First Kind. 

The formal solution of Poisson's equation with source p is given by 

„ 2(F ) = _G/0!^ . (4.6) 
J r - r 



Moreover, the Green's function in toroidal coordinates has the expansion [1 I 

— = N(a,ilj) N(a ,ij ) C "in cos m(cj) - O ) cos n(ip - ip ) 

m,n=0 

xP™ i (cha < )Q^ i (cha > ) , (4.7) 



I r - r o 
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where cr> = max{cr, a }, <j < = min{cx, cr }, and the functions P™_\ and Q m _i are toroidal 
Legendre functions ||12|| . The numerical coefficients are 



Cn 



T(n — m + i 



(4- 



1 {n + m + 

where e n = 2 — <5°. 

In toroidal coordinates, we have 

Q(r -a)5(z) = a' 1 N 2 (a^)5(^) . (4.9) 

Inserting ( |4.9| ) , together with (^4.7|) and ( }4.8|) , into ( |4.6|) , we find the following expression for 
i>2 in toroidal coordinates: 



00 

v 2 = -2GkN(a,ij) £ [P n (cx) P n _i(cha) + G n (<r) Q n _i(ch<r)] e n cosm/> 



where 



n=0 



GJa) 



| rfriV- 1 (r,0)Q n _i(chr) 

(J 

JdrN-\T,0)P n _ l2 ( C hr) 



(4.10) 

(4.11) 
(4.12) 



and we use the abbreviations P„ 1 = P° 1 , Q„ 1 = Q° 1 for the toroidal harmonic 
functions. 



We now turn to the solution of the 77 equation. Writing ( |3.15| ) in full 

8nGk5(r - a) 5(z) , 



d 2 rj d 2 7] 



(4.13) 



dr 2 dz 2 

we can interpret it (after formally replacing r by x) as the potential produced by a uniform 
thin rod lying in the x-y plane, parallel to the 7/-axis at x = a. This is a standard problem 
in electrostatics |13|] , and in cylindrical coordinates, the solution reads 

r](r,z) = 2Gk\og{[(r-a) 2 + z 2 ]/R 2 } , (4.14) 

where Ro is a constant to be discussed later. 

The last remaining metric function to be found is (. Solving ( |3.16| ) with the above 
77, we obtain in cylindrical coordinates: 

2Gk 



C(r,z) 



-Vo + 



log( 



[r - a) 2 + z 2 

Rq 



log( 



[r — a) 2 + z 2 



/ _! ,r — a N 
2r + 2z tan V ) + tan 

v z ' z 



a 2 + z 2 
a , 



(4.15) 
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in which an arbitrary integration function of z has been determined by the requirement that 
the solution be regular on the z-axis. Calculating the angular deficit about the z-axis yields 



5(p = 2vr(l-e~' 70 ) . (4.16) 
Consequently, demanding that there be no conical singularity along the z-axis imposes 

m = o. 

Having determined v, 77, and (, we have formal solutions of the weak-field Einstein 
equations. We must now check whether these solutions are valid within the weak-field 
approximation. It is clear from (|3.14|) , ( |3.15|) that v and 77 are determined only up to 
an additive constant. Moreover, by virtue of (|3.16|) , any constant added to £ may be 



absorbed in a redefinition of r/o, leaving the structure of the equation unchanged. This 
freedom of an additive constant in each of u, rj, and ( enables one to ensure that the weak- 
field approximation is valid at least near the string, the region of spacetime that will be 
relevant when investigating conical singularities and the corresponding deficit angle in the 
next section. (In the standard weak-field approximation in Cartesian coordinates around 
Minkowski space, these constants are naturally chosen so that the potentials vanish at 
infinity.) Furthermore, all future results will be seen to be independent of the particular 
choice of these additive constants. 

Within the model of an infinitely thin string, there is no natural way to determine 
the additive constants, since any choice leads to some of the potentials becoming infinite 
near the string core, invalidating the weak-field approximation. (For instance, irrespective 
of the value of Rq in ([4.14|), 77 is infinite at the ring.) However, a physical string always 



has a non-vanishing thickness, something, strictly speaking, beyond the realm of our thin- 
string model. To give, nevertheless, a value to these constants, we may formally consider 
a slightly fattened string of core radius -R cor e, and fix the additive constants by demanding 
that the potentials be small at distances of order -R core from the infinitely thin ring. (For 
instance, taking Rq = R core in ( |4.14| ) leads to 77 ~ at distances of the order of -R cor e from 



the infinitely thin string.) In what follows, we stay within the infinitely thin string model, 
and we implement the weak- field approximation by the requirement that v, 77, £ << 1 on 
distance scales of order -R cor e from the infinitely thin string. The actual value selected for 
-Rcore does not follow from the thin-string model, and must be supplied by the underlying 
physical theory. We emphasize that we are mathematically dealing with an infinitely thin 
string core, with only the faint memory of its physical origin in the choice of R CO re- It is, 
however, important to stress that all future results will be independent of the particular 
choice of R cnTP - 
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5. Angular Deficit 



To check the metric (|4.2|) for conical singularities on the ring a = oo, we must examine 
circles a = <jq at constant t and <ft around the ring, and calculate the ratio of the proper 
perimeter to the proper radius in the limit that the proper radius tends to zero. For t and 
(f) constant, the metric (|4.2| ) becomes 



ds 2 = a 2 N- A e 2 ^- 2u {da 2 + d^ 2 ) , (5.1) 
and the angular deficit 5ip is given by 



Sip = 2ix — lim 



d^(N- 2 e^)U=a daN~ 2 e^ v 



(5.2) 



From this formula, one sees that the value of dip is independent of any additive constants 
that might appear in v or rj, as emphasised above. In this calculation we are thus free to 
ignore such additive constants, and this will be done below. 



In order to evaluate the limit appearing in (|5.2| ), we require the asymptotic form of 
v and rj for large values of a. The contribution of v\ from (|4.5| ) to the asymptotic behaviour 
of v for a tending to infinity is readily found to be 

ui — > — 2G (/i + k) a , o — > oo . (5.3) 

(In establishing this formula, we use the asymptotic expression K(th(cr/2)) — > a/2 when 
cr tends to infinity JT3| .) In the appendix, we show that v 2 — > as a — >• oo. Thus, 
the asymptotic value of v does not depend on z/ 2 , and this shows that the radial stresses 
T r r = (k/r) 0(r — a) S(z) do not contribute to the potential near the ring. Furthermore, 
the asymptotic behaviour of rj near the ring is 

r}(a,ifj) — > —4Gka as a — > oo . (5.4) 
With flOf) and (|5.4|) , it is straightforward to evaluate the integral (|5.2|) for the deficit 



angle 5^ as 

5^ = 47rG(/i-A;) , (5.5) 
which will be analysed in the conclusion. 

The same result for the angular deficit can be obtained in a different manner, which 
also sheds light on the spacetime structure near the ring in the weak-field limit. Using the 
asymptotic forms of v and rj near the ring, the metric ( [5. 1| ) becomes 

ds 2 = 4a 2 e- 2 ° e +AG ^~ k ^ +2b (da 2 + d^ 2 ) , a -> oo , (5.6) 
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where b denotes the dimensionless combination of all additive constants appearing in the 
potentials. If we now seek coordinates (a', if)'), so that, near the ring, the metric is locally 
flat, 

ds 2 = Aa 2 e" 2a ' (da' 2 + difj' 2 ) , a -> oo , (5.7) 
we see that we must take 

e b e- {{1 - x)a} da = e^' da' , (5.8) 
where A = 2G(/j, — k). Thus, to first order in A, this yields 

a ' = (I- \) a - \-b . (5.9) 

Consequently, if we define ip' — (1 — X)ip, we note that we do obtain the metric in the form 
( |5.7| ). It is then clear that if)' runs from to 2ir(l — A), and that we thus recover the deficit 
angle 

5if) = 4nG(n-k) (5.10) 

already found in (|5.5| ), once again seen to be independent of the choice of the additive 
constants. A similar method was employed by Vilenkin [[|. 



6. Conclusion 

In this paper, we constructed a stationary, axially symmetric metric satisfying Einstein's 
weak-field equations for a source describing an infinitely thin ring of radius a, linear mass 
density fi, and azimuthal stress per unit length k. The conservation of stress-energy required 
the presence of a radial stress 

T\ = -e(r-a)6(z) (6.1) 
r 

to support the ring. The metric obtained exhibits a conical singularity along the ring, with 
deficit angle 

dip = AnG(fi-k) . (6.2) 

In the case of a ring of pressureless dust (k — 0), a conical singularity arises, with 
deficit angle 4nGfi. This is usually interpreted by saying that, because of the presence of 
the singularity, such a source cannot remain stationary and must collapse. 

When the ring is made of "ordinary" matter (k > , k << /i), the conical 
singularity is still present but is less severe than for pressureless dust. Such a source is also 
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considered as non-physical since, unless the (positive) pressure equals the energy density, 
the angular deficit 8ip is still non-zero. 

When the ring is made of "string" matter (k = —fi), we recover Vilenkin's result 
5ip = 8nGfi. This result is interesting for two reasons: first, it shows that one half of 
the angular deficit of the string source is of "non-string " origin (since one half of the effect 
remains for pressureless dust). Secondly, it establishes from the field equations, within the 
thin-string model, that indeed, a conical singularity is present along a circular string, and 
that the angular deficit takes the same value as for a straight string of the same linear 
mass density. This was an assumption made by FIU in their investigation of circular cosmic 
strings at an instant of time symmetry. Our formalism provides a weak-field proof of the 
validity of the FIU hypothesis. 
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Appendix A. Asymptotic Behaviour of z/2 near the Ring 



Establishing the asymptotic behaviour of the potential (|4.10| ) near the ring, namely, for 



a tending to infinity, is a delicate matter. This is due to the fact that the source, being of 
the form (k/r) @(r — a) S(z), decreases very slowly at infinity when the radius r increases 
without limit. The difficulty at infinity in r translates into a difficulty at in a. 

For the calculation of F n (a) of ( |4.11| ), however, the point o = does not belong to 
the domain of integration and one finds, after replacing Q n _i(cho-) by its asymptotic value 
T5| and performing the integration: 



(c ) - (2n) 1 / 2 T(n + ^)e^ n+1 ^/(n + l)\ . (A.1) 



On the other hand, a problem does arise in the calculation of G n (a): Near r = 0, the toroidal 
harmonics P n _i (chr) tend to 1 for all n | I6|| , whereas the factor iV _1 (T, 0) appearing in the 
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integrand of ( |4.12 ) tends to 2 1 / 2 r~ 1 . As a result, G n diverges, and a cut-off e must be 
introduced. If this is done, the integrand is finite everywhere, and the regularised G n (a), 
denoted by G^(a), reads: 

(A.2) 



G» ee rfriV- 1 (r,0)P n _i(chr) . 

To go further, one must consider the asymptotic value of P„_i(chcr). 

For n — 0, the corresponding toroidal harmonic is expressible in terms of the Com- 
plete Elliptic Integral of the First Kind [117 



as 



1 



1 



P_i(cha) = 27r- 1 K(th(-a))/ch(-a) . 

For a very large, this is asymptotically given [|14]] by 

P_i(ch(j) -> 2n- 1 ae- a/2 . 
For n > 1, the asymptotic expression []TH| of P n _i(chcr) is 

P n _i -> 7r- 1/2 (n-l)!e {n - 1/2) 7r(n + i) , n>l 

2 



(A.3) 



(A.4) 



(A.5) 



When the asymptotic expressions ( |A.4| ), ( |A.5| ) are inserted into the integral ( A.2 ) for 



G„ (er), the following results are obtained: 



G » 



G?(a) 



A(e) , A(e) 



dr N~ (t } 0)P_i(chr) 



-) 1/2 -/r(|) 

7T 2 

-) 1/2 (n-2)!e (n - 1) 7r(n+ J) , n>2 

7T 2 



so that, finally, the asymptotic form of the series fl4.10|) becomes 



V-2 



-Gk 



2 1 / 2 7rA(e) + o(e- a ) 



(A.6) 
(A.7) 
(A.8) 

(A.9) 



on 



As one can see, for every non- vanishing value of e, the potential u 2 tends to a constant 
the ring, that is, for a tending to infinity. The fact that v 2 is asymptotically constant is 
sufficient to establish the result ( |5.5| ) of the deficit angle, irrespectively of the value of the 
constant. It is, however, possible to go one step further and renormalise the potential so 
that z/ 2 vanishes at the ring. This amounts to subtracting from z/ 2 the (a, ^-independent 
term 2 l l 2 i\A(e). (Although, when the cut-off e tends to 0, the constant A(e) diverges, this 
procedure is nothing more than a conventional renormalisation.) 
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